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f  the  flow  through  the  passage  and  includes  a  nonlinear  description  of  cascade  loss  anc 
turning  correlations  from  available  experimental  sources.  The  flow  field  analysis 
ir,vOi''„*s  the  nonlinear,  time-dependent  equations  for  the  vorticity  and  the  stream 
function.  Coupling  of  the  two  is  accomplished  through  the  boundary  conditions  by 
r.utual  relationships  between  the  pressure  change  across  the  blade  row  and  the  change 
in  vorticity  in  the  flow  field  analyses. 

Within  the  assumptions  that  the  flow  is  two-dimensional  and  incompressible,  the 
umcrical  solution  is  capable  of  predicting  the  influence  of  an  upstream  axial  distor¬ 
tion  on  the  onset  of  a  circumferentially  rotating  stall  pattern  for  a  single  blade  row. 
The  speed  of  rotation  of  the  stall  cell  and  the  spatial  attenuation  of  the  distortion 
wave  are  also  predicted,  and  although  the  observed  experimental  data  are  generally 
available  only  for  multistage  systems,  the  predicted  results  for  the  single  blade  row 
are  in  qualitively  good  agreement  with  the  data. 

The  two  major  parts  of  the  theory  can  be  decoupled  and  the  blade  row  analysis  can 
uc  performed  separately.  For  a  prescribed  axial  distortion  measured  at  the  blade  row 
inlet  face,  this  portion  of  the  theory  predicts  a  circumferential  pressure  distribution 
it  the  blade  row  exit  which  is  in  good  agreement  with  measured  exit  data. 
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Unsteady  Fluid  Dynamic  Response  of  an 
Axial-Flow  Compressor  Stage  with  Distorted  Inflow 


SUMMARY 


A  nonlinear,  large  disturbance  theory  has  been  developed  which  couples, 
interactively,  the  flow  through  the  blade  passages  of  a  turbcsiachine  blade 
row  and  an  axially  distorted  flow  field.  The  blade  row  analysis  is  based 
on  the  time-dependent  energy  equation  of  the  flow  through  the  passage  and 
includes  a  nonlinear  description  of  cascade  loss  and  turning  correlations 
from  available  experimental  sources.  The  flow  field  analysis  involves  the 
nonlinear,  time-dependent  equations  for  the  vorticity  and  the  stream  function. 
Coupling  of  the  two  is  accomplished  through  the  boundary  conditions  by 
mutual  relationships  between  the  pressure  change  across  the  blade  row  and 
the  change  in  vorticity  in  the  flow  field  analyses. 

Within  the  assumptions  that  the  flow  is  two-dimensional  and  incompress¬ 
ible,  the  numerical  solution  is  capable  of  predicting  the  influence  of  an 
upstream  axiol  distortion  on  the  onset  of  a  circumferentially  rotating  stall 
pattern  for  a  single  blade  row.  The  speed  of  rotation  of  the  stall  cell  and 
the  spatial  attenuation  of  the  distortion  wave  are  also  predicted,  and 
although  the  observed  experimental  data  are  generally  available  only  for 
multistage  systems,  the  predicted  results  for  the  single  blade  row  are  in 
qualitively  good  agreement  with  the  data. 
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The  two  major  parte  of  the  theory  can  "be  decoupled  and  the  "blade  row 


analysis  can  be  performed  separately.  For  a  prescribed  axial  distortion 
measured  at  the  blade  row  inlet  face,  this  portion  of  the  theory  predicts 
a  circumferential  pressure  distribution  at  the  blade  row  exit  which  is  in 
good  agreement  with  measured  exit  data. 


INTRODUCTION 


At  Its  inception  the  object  of  this  program  was  to  assess  the  extent  to 

which  the  performance  characteristics  of  a  compressor  operating  in  a  circum- 
0 

ferentially  distorted  inlet  flow  are  influenced  by  unsteady  aerodynamic 
effects.  A  first  effort  (Ref.  l)  concentrated  on  the  development  of  an 
unsteady  model  based  on  unsteady  aerodynamic  data  from  an  isolated  airfoil, 
and  it  was  shown  that  a  reasonable  correlation  could  be  attained  between 
measured  and  predicted  circumferential  pressure  variations.  Although  this 
model  yielded  encouraging  results  it  was  recognized  that  it  depended  strongly 
on  a  number  of  intermediate  transformations  and  assumptions  to  convert 
isolated  airfoil  data  to  a  form  compatible  with  compressor  operation.  A 
more  serious  limitation  of  the  model  in  Ref.  1  involved  the  use  of  measured 
upstream  flow  properties  obtained  at  the  inlet  plane  during  an  inlet  distor¬ 
tion  experiment.  These  upstream  conditions  are  strongly  dependent  on  the 
presence  of  the  blade  row,  as  shown  in  Fig.  1.  Here  the  measured  stall  angle 
parameter,  a  ,  has  been  plotted  versus  circumferential  position,  8  ,  for  a 
180  deg  distortion  condition  analyzed  in  Ref.  1.  The  dashed  line  in  Fig.  1 
represents  the  angle  of  attack  distribution  that  would  be  measured  immedi¬ 
ately  behind  the  distortion  screen  far  upstream  of  the  blade  row,  and  the 
solid  line  is  the  actual  angle  of  attack  distribution  at  the  blade  row  inlet. 
Is  is  seen  that  the  mutual  interaction  between  the  blade  row  and  the  circum¬ 
ferential  velocity  profile  far  upstream  of  the  blade  row  produces  a  further 
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distortion  of  large  magnitude  on  the  previously  distorted  flow.  However, 
at  the  conclusion  of  the  work  reported  in  Ref.  1  it  was  not  possible  to 
predict  the  change  in  the  distortion  pattern  as  it  came  under  the  influence 
of  the  blade  row,  and  it  was  necessary  to  measure  the  upstream  conditions 
during  an  actual  distortion  experiment. 

One  of  the  principal  objectives  of  the  present  work  is  to  provide  the 
engine  designer  with  the  analytical  means  to  estimate  the  distortion  toler¬ 
ance  or  the  distorted  flow  performance  of  a  rotor,  including  unsteady  effects, 
without  resorting  to  an  expensive  test  program.  This  means  the  designer  must 
be  able  to  prescribe  a  distortion  far  upstream  of  the  rotor  face,  and  the 
analysis  must  be  capable  of  properly  modifying  the  distortion  at  the  rotor 
race  in  accordance  with  the  coupling  between  the  upstream  flow  distortion 
and  the  downstream  boundary  conditions. 

It  was  just  shown  (cf.  Pig.  l)  that  the  coupling  between  the  downstream 
flow  and  the  upstream  distortion  is  quite  strong,  and  it  is  not  possible  to 
prescribe  the  conditions  at  the  compressor  face  a  priori  or>  the  basis  of 
upstream  distortion  profile  and  geometry  alone.  Hence  an  effort  was  made  to 
derive  an  analytical  method  to  predict  the  behavior  of  the  entire  coupled 
flow  field.  Several  analyses  (Refs.  2,  i,  4)  hav?  attempted  to  do  this 
using  linearized  theory  (based  on  the  assumption  of  a  small  disturbance 
compared  to  the  mean  velocity).  However,  this  restriction  is  inappropriate 
for  the  inlet  distortion  problem  which  hao  large  amplitude  flow  fluctuations 
which,  could  cause  the  compressor  to  operate  beyond  the  steady-state  stall 
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boundary  over  a  olgniflrnnh  cJ rcumferiHitial  orctor.  Accord  Ingly ,  Uie 
present  method  involves  a  numerical  solution  of  a  set  of  nonlinear  equa¬ 
tions  which  are  not  restricted  by  a  small  disturbance  requirement,  and  which 
account  for  the  interactive  influences  of  the  blade  row  and  the  entire 
d istorted  flow  field . 

In  particular,  the  present  stability  analysis  was  initially  divided 
into  two  parts  which  modelled  the  blade  row  and  the  flow  field  separately. 

Uy  this  means  it  was  possible  to  demonstrate  the  validity  (and  tractability ) 
of  each  constituent  part  before  attempting  to  couple  them  in  a  single 
stability  analysis. 

Tn  the  first  part  the  earlier  blade  row  model  of  Ref.  1  was  replaced 
with  a  more  direct  representation  which  accounted  for  the  presence  of 
neighboring  blades.  This  new  model  describes  the  unsteady  flow  through  a 
biade  passage  of  either  a  rotor  or  a  stator  blade  row.  It  is  based  on  the 
time-dependent  energy  equation  of  the  flow  through  the  passage  and  includes 
a  nonlinear  description  of  cascade  loss  and  turning  correlations  from 
available  experimental  sources.  Use  of  this  model  with  prescribed  upstream 
''•orid itionn  at  the  blade  row  inlet  plane  yielded  good  agreement  with  measured 
data . 

The  second  part  of  ohe  present  analysis  was  concerned  with  developing 
a  flow  field  model  for  the  regions  upstream  and  downstream  of  the  blade  row. 
This  model  was  designed  to  be  capable  of  accounting  for  mutual  interactions 


between  the  blade  row  and  the  upstream  inlet  distortion.  This  part  consists 
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of  a  nonlinear,  two-dimensional,  time  dependent  flow  field  analysis 
involving  the  equations  for  the  vorticlty  and  the  stream  function. 

Initially  this  flow  field  solution  was  decoupled  from  the  blade  row 
analysis  to  provide  a  baseline  transient  response  problem  suitable  for 
checking  out  the  stability  of  the  required  numerical  procedures. 

Finally,  the  two  parts  of  the  analysis  were  coupled.  The  upstream 
8nd  downstream  flow  transients  contributed  to  the  time-dependent  blade 
passage  flows,  and  the  unsteady  blade  row  pressure  rise  provided  a  mid¬ 
flow  boundary  condition  for  the  time-dependent  flow  field  analysis.  Direct 
comparisons  with  experiment  were  not  made  because  of  the  lack  of  experimental 
data  on  the  response  of  an  isolated  rotor  subjected  to  an  upstream  distor¬ 
tion.  Nevertheless,  it  is  felt  that  the  results  show  reasonable  trends 
which  are  qualitatively  correct  when  compared  with  available  experimental 
data  for  a  multi-stage  compressor. 
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BLADE  ROW  ANALYSIS 

Consider  the  flow  of  an  incompressible  fluid  through  one  of  the  blade 
passages  shown  in  Pig.  2.  The  flow  is  assumed  to  be  two-dimensional,  and 
it  enters  the  passage  at  station  1  and  exits  at  station  2.  Continuity 
requires  that  the  flow  which  exists  at  station  2  must  be  equal  to  the  flow 
entering  at  station  1.  This  can  be  expressed  by  the  equation 


fydU,dy  1  U.t  = /y° U2dy  =  U2r  £Q 

\  \ 


(1) 


where  Q  is  the  flow  rate  through  the  passage,  and  the  bar  over  the  variables 
Ui  and  U2  denotes  gap-averaged  quantities.  Thus,  the  continuity  equation 
requires  that  the  average  axial  velocity  be  a  constant  at  every  point  within 
the  blade  passage 


u,  =  u2--u 


(2) 


In  the  present  problem  the  flow  through  each  blade  passage  is,  in 
general,  unsteady.  Hence  the  accompanying  variation  in  energy  flux  through 
each  passage  will  be  accounted  for  using  the  unsteady  energy  equation  for  an 
incompressible  fluid  (cf.  Eq.  (2.2l)  in  Ref.  5)  flowing  through  a  control 
volume-  When  this  equation  is  rewritten  in  the  present  notation  we  obtain 

/  /r  (?/,q2)dv*/  (p  +-h2H*ncJA  =  0  (3) 

v  ate  do  ' 


(3) 


where  v  in  the  first  integral  denotes  the  control  volume  and  the  notation 
abcda  under  the  line  integral  denotes  the  control  surface  in  Pig.  2  over 
which  it  is  to  be  evaluated.  The  line  integrals  along  the  airfoils,  ab  and 
cd,  are  zero  because  the  flow  is  everywhere  tangent  to  these  surfaces,  so 
q.ndA  =  0.  The  remaining  integrals  over  be  and  da  may  be  simplified  by  using 


q,*ndA  :  -U,dy  ) 
q2*ndAxu2dy  ) 


so  Eq.  (l)  becomes 


fyc(p2+  T/>qf)u2dy -/*!>,  +  ypq,2)u,dy  +  /  ~  (^/>q2)dv=o 

yb  Vq  v 


Bernoulli's  equation  may  now  be  used  to  replace  the  sum  of  the  static  and 
dynamic  pressures  at  the  inlet  and  exit  planes  by  total  pressures. 


fC  pQ2My-/0P0,My  =  - 1  af  (i^2)dv 


The  validity  of  the  next  portion  of  this  analysis  depends  on  our  ability 
to  perform  meaningful  gap  averages  of  many  of  the  variables.  Inherent  in 
thi6  requirement  is  the  concurrent  requirement  that  these  variables  be,  in 
general,  slowly  varying  circumferentially  with  wavelengths  that  are  greater 
than  blade  gap.  Under  these  conditions  the  variable  will  simply  be  replaced 
by  its  local  value.  For  example,  if  both  the  inlet  total  pressure  and  the 
inlet  velocity  meet  this  requirement,  the  second  integral  in  Eq.  (6)  may  be 
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rewritten  ao 


jfr  pOi  uidy  =  po,Ur  =  PolQ 


(?) 


If  the  assumption  of  no  losses  in  the  flow  were  made,  the  same  result  would 
he  obtained  for  the  first  integral  in  Eq.  (6).  Instead,  it  will  be  assumed 
that  a  total  pressure  loss  occurs  that  is  associated  with  blade  wakes  and 
Pq2  and  1^2  will  be  replaced  by  P02  +  AP02  and  Xfe  +  Alfe,  respectively.  Here 
the  bar  represents  the  gap  average  of  a  quantity  which  varies  rapidly  in  a 
local  region,  such  as  the  wake,  but  which  has  an  otherwise  slowly  varying 
characteristic.  The  variables  APQg  and  are  local  deviations  from  the 
gap-averaged  quantities  and  represent  some  measure  of  the  loss,  to  be 
resolved  in  the  analysis  below.  The  first  integral  of  Eq.  (6)  will  now  be 
replaced  by 

/  p02u2^y  z  f  ^P02  +  <ApOz)  (U2+  AU2)dy  =  Pq2U2  T  +  f  CAPQ2AU2dy 


~  po202r  +  TPQ,2X  U2r 


(8) 


=  *02^  *  ?  Pd 1^  X  Q 


where  X  is  an  unsteady  loss  parameter  (see  below). 

The  unsteady  loss  parameter,  X,  is  conceptually  defined  in  Eqs .  (8)  as 
the  instantaneous  gap-averaged  deviation  of  exit  total  pressure  during  an 
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unsteady  process.  In  a  steady-state  proces-  a  combination  Oi  Eqs.  (6), 
(T),  and  (8)  shows  that  X  reverts  to  its  traditional  steady-state 
definition, 


P0l"P02 

i/>pi2 


(9) 


Although  the  specific  unsteady  behavior  of  X  ia  not  known,  experience  has 
shown  that  unsteady  processes  do  not  follow  the  steady  characteristics,  but 
Jag  them  by  some  appropriate  time  constant.  Consequently,  a  hypothetical 
first  order  lag  model  will  be  chosen,  as  in  Refs.  6  and  7> 


coxaT 

U* 


dX 

dt  ' 


+  X 


(10) 


where  a ^  is  a  dimensionless  time  constant  associated  with  the  chordwise 
convection  of  disturbances  within  the  boundary  layer  and  where  U*  is  the 
circumferential  average  of  the  axial  velocity. 


U 


»  _ 


(11) 


(Note  that  8  and  y  are  both  tangential  coordinates.)  This  hypothetical 
model  will  be  used  presently. 

Consider  next  the  rate  of  change  of  kinetic  energy  within  the  passage, 
represented  by  the  right  hand  side  of  Eq.  (6).  The  square  of  the  local 
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velocity  can  be  rewritten  sb 


q2  =  U2+V2  =  U "( 1  +  ton2£)  =  U2sec2C 


(12) 


where  §(x,y,t)  is  the  local  flow  angle  within  the  passage.  Hence,  the 
volume  integral  becomes 


f 


_d 

at 


2  pwl  L.,  usec£d''dx 


YU) 


(13) 


where  y^(x)  and  y2(x)  are  the  y-coordinates  of  blades  ab  and  dc,  respectively 
(cf.  Fig.  2).  The  evaluation  of  this  integral  depends  on  a  number  of 
assumptions.  First  it  is  assumed  that  the  flow  within  each  passage  is  mostly 
one -dimensional.  In  other  words,  the  assumption  is  made  that  variations 
normal  to  the  mean  streamline  through  the  passage  are  small,  which  implies 
that  separate  gap  averages  of  the  integrand  factors  can  be  taken,  or 

rc ox  /*yaW  g  2+  rcax  ~z  £ 

J  I  u  sec  C  dy  dx  =  rj  u  sec  £  dx  (14 ) 

*0  y,(x)  •'o 


From  wvove,  %  is  independent  of  x  through  the  blade  passage.  Furthermore, 
it  will  be  assumed  that  U  is  slowly  varying  in  the  circumferential  direction, 
and  in  accordance  with  the  convention  established  earlier,  the  gap-averaged 


&j&’i&&KAMt m+m  »..  - 


removed  from  under  the  x-integral  in  Eq.  (13),  which  now  becomes 


{  ir  (ip?) dv  =  jpt  af  °T0<  sec2  £  d* 


Finally,  the  assumption  of  small  gap  to  chord  ratio  will  be  made.  Under 
this  assumption  virtually  all  of  the  flow  angle  changes  will  take  place  in 
the  entrance  region  of  each  passage,  the  chordvise  extent  of  which  is  small 
compared  with  the  chord,  and  the  flow  through  the  remainder  of  each  passage 
will  behave  like  a  channel  flow  (Ref.  8).  Hence,  the  chordwise  integral  of 
sec5? in  Eq.  (15 )  can  be  replaced  by 


jfCoxsec2£dx  s  jQCQ*sec2£me,oldx  5  c0~m 


where  5  metal  ls  the  angle  of  the  blade  camber  line,  commonly  called  the 
metal  angle,  and  where  r-m  now  represents  a  chordwise  average.  Furthermore, 
this  quantity  is  independent  of  time  by  virtue  of  its  channel-like  behavior, 
so  Eq.  (15)  becomes 


l  IT  (i'>q2)du  =  ?Pr  CoxSmx  OJ =PQc0IEm|y- 


When  Eas .  (7).  (8),  and  (17)  are  substituted  into  Eq.  (6)  the  result  is 


p02~  Pot  =  AP0=  “  jp<\?X-p  cox«m~^- 
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(which  is  equivalent  to  Eq.  (9)  in  Kef.  7)*  Once  again  Bernoulli's  equation 
will  be  used  at  both  inlet  and  axit  planes, 

poi=  Pi  +  ^PQ  i2 

P02  =  P2  +  j[P  $2 

Here  a  specific  notation  for  gap  averaging  is  not  necessary  on  P2  because  in 
the  definition  of  P02  (Eas.  (8))  all  of  the  local  wake  deficiencies  have  been 

lumped  together  in  the  loss  parameter,  X.  Hence,  both  P2  an^  q2  can  be 
regarded  as  slowly  varying  quantities  and  can  be  represented  by  local  values. 
It  was  stated  earlier  that  the  equation  of  continuity  requires  that  the  gap- 
averaged  axial  velocity,  U,  be  constant  through  the  blade  row.  Thus  we  may 
write 

U  =  q,  sin  =  q2sm/32  (20) 


Again,  no  specific  gap  averaging  notation  is  necessary  for  02  because  its 
value  will  be  taken  from  a  set  of  gap-averaged  correlation  results  such  as 
those  in  Ref.  9*  With  the  introduction  of  Eqs .  (19)  and  (20),  Eq.  (18)  now 
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APPLICATION  OF  BLADE  ROW  ANALYSIS 


Before  the  blade  row  analysis  was  coupled  with  the  flow  field  analysis 
(see  subsequent  sections  below)  it  was  decided  to  test  its  validity  by  using 
it  to  examine  some  of  the  experimental  data  previously  analysed  in  Ref.  1. 

In  this  application  the  measured  inlet  conditions  at  the  rotor  face  were 
prescribed  to  determine  rotor  exit  conditions ,  end  the  measured  inlet  condi¬ 
tions  at  the  stator  face  (i.e.,  the  rotor  exit  conditions}  were  prescribed 
to  determine  stator  exit  conditions.  All  of  the  quantities  on  the  right  hand 
side  of  Eq.  (21 )  were  known  or  easily  derived.  For  example,  steady-state 
loss  and  turning  characteristics  for  the  rotor  and  stator  had  been  previously 
determined  experimentally  and  are  shown  in  Fig.  3  versus  inlet  angle,  Bi« 

The  time  derivative  in  Eq.  (21 )  was  rewritten  in  terms  of  rotor  rotation 
speed,  d/dt  =  0d/d9  ,  and  the  loss  coefficient  time  constant  in  Eq.  (10) 
was  assigned  the  value  orqr  =  1,  in  accordance  with  the  measurements  of  Ref. 
10.  The  inlet  angle,  ,  was  then  determined  as  a  function  of  circumferen¬ 
tial  position,  6  ,  and  Eq.  (10)  was  solved  for  X(  0  )  using  Xss  (  R1(  8  ))  as  a 
known  input.  Finally,  the  exit  static  pressures  were  calculated  for  both 
rotor  and  stator  using  Eq.  (21 ). 

Three  specific  cases  for  l80  deg  circumferential  screens  from  Ref.  1 
were  examined  and  are  presented  in  Figs,  k,  and  6  for  the  mid  flow,  peak 
pressure,  and  surgp  conditions,  respectively  (cf.  Figs.  18,  9,  and  22  of 
Ref.  l).  In  each  these  figures  the  square,  triangular,  and  circular 
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symbols  denote  measured  rotor  inlet,  rotor  exit,  and  stator  exit  pressures. 

The  solid  curves  arc  the  exit  pressures  predicted  by  Eq.  (2l).  A  comparison 
with  the  results  in  Ref.  1  showed  that  Eq .  (21)  generally  improved  the  cor¬ 
relation  with  the  experimental  data.  Ihi6  is  particularly  true  in  the 
prediction  of  the  rotor  exit  pressure  for  the  peak  pressure  condition  as 
shown  in  Fig.  7« 

The  apparent  discrepancies  for  the  stator  prediction  in  all  three  cases 
can  be  ascribed  to  a  general  inadequacy  of  the  steady-state  stator  performance 
data  which  had  to  be  extrapolated  to  encompass  the  unsteady  operating  limits 
during  distorted  operation.  The  disagreement  between  predicted  and  measured 
rotor  pressure  in  Fig.  6  at  the  screen  exit  (  0  »  90  deg)  is  also  believed 
to  be  ascribable  in  part  to  a  required  extrapolation  of  r'  performance 
data.  More  importantly,  the  combination  of  deep  stall  penetration  and 
unsteady  flow  at  the  screen  edge  probr.bly  caused  a  flow  instability  that  was 
not  measureable  it:  detail  with  +h«*  steady-state  instrumentation  used  in  this 
experiment.  Furthermore,  use  of  Eq.  (21 )  with  spatially  fixed  input  variables 
can  only  yield  spatially  fixed  exit  predictions,  and  it  will  be  shown  in  a 
later  section  that  a  coupling  between  the  flow  field  equations  and  Eq.  (21 ) 
leads  to  both  temporal  and  spatial  flow  instabilities  unasr  conditions  of 
deep  stall  penetration.  Crusequently  it  was  felt  that  most  of  the  observed 
discrepancies  in  Figs.  4,  and  6  could  be  discounted,  and  that  the  agreement 
between  rotor  exiv  pressure  prediction  and  experiment,  shown  in  Fig.  5>  was 
significant  and  provided  the  required  validation  of  Eq.  (21 ). 


A  closer  examination  of  the  effects  of  individual  terms  in  Eq.  (21 ) 
was  performed  and  the  results  are  shovm  in  Fig.  8.  Ab  before  the  measured 
data  are  denoted  by  the  circled  points,  and  the  theory,  using  Eq.  (21 )  with 
Qfrp  =  1,  is  shown  by  the  solid  line.  The  quasi-steady  theory,  using  the 
cascade  correlations  of  Fig.  3>  is  shown  as  a  dashed  line.  This  prediction 
was  obtained  from  Eq.  (2l)  by  setting  X  =  Xss  and  pCax  JTmdU/  St  m  0. 
Finally,  a  theoretical  curve  including  unsteady  loss  but  neglecting  the 
kinetic  energy  of  the  flow  within  the  passage,  denoted  by  the  dash-dot 
curve,  was  obtained  by  setting  =  1  in  Eq.  (10)  and  suppressing  the  last 
term  of  Eq.  (21 ).  Figure  8  shows  that  unsteady  terms  must  be  included  in  the 
anrlysis  to  properly  model  the  fluid  dynamics  of  inlet  distortion.  It  also 
ohows  that  for  the  specific  case  examined  here  the  unsteady  kinetic  energy 
term  has  a  greater  influence  on  the  results  than  the  unsteady  loss.  (The 
difference  between  the  dashed  and  the  dash-dot  curves  is  the  contribution  of 
the  unsteady  loss,  and  the  difference  between  the  dash-dot  and  the  solid 
curves  is  that  of  the  unsteady  kinetic  energy  within  the  blade  passage . ) 
However,  it  will  be  shown  in  a  later  section  that  the  unsteady  loss  term  has 
a  significant  effect  on  flow  stability  and  rotating  stall. 
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behind  an  axial  compressor  blade  row,  unwrapped  tangentially  and  projected  \ 

I 

onto  a  two-dimensional  plane.  The  lines  ABC  and  DEF  represent  a  common 
axial  trace  in  the  three-dimensional  compressor  system,  and  this  correspon¬ 
dence  imposes  a  periodic  boundary  condition  on  the  flow;  viz.,  it  is  required 
t.hat  any  flow  particle  leaving  (or  entering)  the  region  through  line  ABC  must 
enter  (or  leave)  the  region  through  line  DEF  with  the  same  vector  velocity. 

A  distorted  inlet  flow  is  prescribed  along  the  inlet  boundary  AD,  and  a 

constant  static  exit  pressure  condition  is  imposed  at  exit  boundary  CF.  The 

blade  row  represented  by  the  line  BE  is  modelled  as  a  semi-actuator  disk 

across  which  there  is  a  change  in  the  flow  direction  as  well  as  a  change  in 

the  total  pressure .  \ 

i 

i 

i 

* 

i 

% 

i 
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The  flow  field  upstream  and  downstream  of  the  blade  row  will  satisfy 
the  incompressible,  time-dependent,  two-dimensional  Euler  equations  of 
motion, 


vorticity  can  only  be  created  at  the  boundaries  of  the  designated  flew 
field.  Thus,  with  reference  to  Fig.  9,  vorticity  can  only  be  created 
along  the  inlet  plane  boundary,  AD,  the  actuator  disk,  BE,  and  the  exit- 
plane  boundary  CF.  In  particular,  the  creation  of  vorticity  associated 
with  reverse  flow  at  the  actuator  disk,  BE,  is  admissible  in  this  solution. 
This  means  that  the  rotating  stall  phenomenon,  which  can  be  accompanied  by 
local  regions  of  reverse  flow,  is  not  excluded  from  thiB  analytical  model. 
The  creation  of  vorticity  at  this  boundary  is  governed  by  the  interactions 
of  the  flow  field  equations  developed  in  this  section  and  the  pressure  rise 
boundary  condition  across  BE  determined  in  an  earlier  section.  It  will  be 
shown  later  that  the  specification  of  the  form  of  the  inlet  distortion  at 
the  entrance  plane  AD  provides  the  necessary  information  for  the  specifica¬ 
tion  of  vorticity  along  this  boundary.  It  is  not  possible  to  establish  any 
rational  cirterion  for  the  specification  of  vorticity  at  the  exit  plane. 

This  would  be  necessary  only  if  a  region  of  reverse  flow  at  the  actuator 
disk  BE  extended  downstream  far  enough  to  penetrate  the  boundary  CF.  Under 
these  circumstances  either  the  case  under  consideration  could  be  regarded  as 
invalid  or  the  boundary  CF  could  be  moved  further  downstream  to  avoid  the 
reverse  flow  region.  However,  in  all  of  the  cases  considered  herein  there 
were  no  detectable  regions  of  reverse  flow  snd  it  is  considered  unlikely  that 
such  problems  will  arise  in  the  vicinity  of  the  exit  plane  boundary. 


The  assumption  that  the  flow  field  is  incompressible  and  two-dimensional 
in  the  cascade  plane  allows  the  introduction  of  a  stream  function,  &  ,  which 
is  related  to  the  axial  and  tangential  velocity  components  by  the  equations 


,V  - 


d\p 

dx 


(27) 


A  combination  of  Eqs.  (25)  and  (27)  yields  the  expression 


t2±  .  , 


(28) 


The  final  flow  field  equations  that  will  be  needed  for  this  solution 
involve  the  relationships  between  the  total  pressure  at  any  point  in  the 
flow  field  and  the  other  fluid  properties  at  that  point.  From  Eqs.  (22) 
through  (25)  (or  from  Ref.  11 )  these  relationships  take  the  form 


<3u  ]  d 
dt  dx 


(29) 


dv  _d_ 
dt  +  dy 


(3C) 


which  are  essentially  the  axial  and  tangential  momentum  equations. 

A  simultaneous  solution  of  Eqs.  (26)  and  (28)  which  satisfies  the 
appropriate  boundary  conditions  will  yield  the  velocity  field  upstream  and 
downstream  of  the  blade  row,  and  the  total  pressure  field  can  be  obtained  by 
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integrating  Eqs.  (29)  and  (30).  When  these  equations  are  combined  with 
Eqs.  (19)  and  (21)  for  the  blade  row  pressure  rise  the  result  comprises 
the  desired  nonlinear,  coupled  set  of  equations,  valid  for  large  disturbances, 
and  including  nonlinear  loss  and  unsteady  flow  effects. 
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NUMERICAL  PROCEDURE  AND  BOUNDARY  CONDITIONS 


The  numerical  procedure  employed  in  this  study  is  a  standard  center 
differencing  method  for  the  stream  function  equation,  (28),  and  Arkawa’s 
vorticity  conservation  algorithm  (Ref.  12)  and  the  Crank-Nicholson  scheme 
(Ref.  13)  to  solve  Eq.  (26).  The  usual  procedure  is  to  prescribe  a  distor¬ 
tion  at  t  *  0  at  the  inlet  plane,  AD,  and  to  maintain  it  at  the  inlet  for 
all  later  times.  This  implies  a  tacit  assumption  that  the  inlet  plane  is 
far  enough  upstream  to  be  unaffected  by  the  presence  of  the  blade  row  at 
BF.  At  each  time  step  over  t  2  0  the  propagation  of  the  disturbance  through 
the  flow  field  is  computed  subject  to  a  simultaneous  satisfaction  of  all 
boundary  conditions .  The  computation  is  terminated  either  when  a  steady 
state  condition  is  achieved  or,  in  the  case  of  an  unsteady  downstream 
solution,  when  such  a  solution  becomes  periodic.  (The  latter  is  characteris¬ 
tic  of  rotating  stall.)  Each  of  the  boundary  conditions  will  now  be  des¬ 
cribed  in  detail. 

Inlet  Conditions 

The  circumferential  distortion  profile  is  prescribed  at  the  inlet  plane 
AD  by  specifying  the  functional  form  of  the  stream  function, 


Note  that  thie  specification  is  one-dimenaional  (i.e.,  independent  of  x  at 
this  point)  and  that  it  can  have  a  time  dependency.  From  Eqe.  (28)  and  (3 l)> 
the  vorticity  along  AT)  can  be  computed  using 


;(x,y,t) 


d2^(x,y,t) 


'xOt  AD 


x  Ot  AO 


=  f(y.t) 


Flow  Bsriodicity 


It  was  noted  in  an  earlier  section  that  flow  entering  (or  leaving) 
along  the  line  ABC  must  disappear  (or  reappear)  in  the  same  form  and  dis¬ 
tribution  along  the  line  DEF.  This  is  expressed  by  the  periodicity  condi¬ 
tions 


cu.y,t)! 


'y  ot  OEF 


ly  Ot  ABC 


'//(x,y,t)  =  ^U,y,t)  to* 

y  Ot  DEF  y  ot  ABC 

where  Q*  is  the  axial  flow  rate  through  the  entire  flow  field, 


w  -  clfn^u 


where  Rpjj  is  the  mean  radius  of  the  annulus . 


Pressure  Rise  Across  Actuator  Disk 


The  expression  for  the  pressure  rise  across  the  actuator  disk,  "based 
on  unsteady  "blade  passage  characteristics,  iB  given  in  Eq.  (l8).  A 
comparable  equation  for  total  pressure  across  the  blade  row  may  be  obtained 
by  rewriting  Eq.  (30)  in  the  form 


When  Eq.  (l8)  is  differentiated  with  respect  to  y  it  can  be  combined  with 
Eq*  (36)  to  eliminate  the  pressures,  and  the  result  is 


’Mr--fr+uM.) 


where  U  »  =  U2  by  continuity.  Analogous  to  Eq.  (20)  (or  from  Fig.  2)  we 

can  write 


so  Eq.  (37)  finally  becomes 


4£-4-/oton  B.) 

C/i  Oi  \  •  * 


“T'4"(q*X)“C«Sn>‘^5t‘  +  U(^  £*) 


This  equation  relates  the  velocity  field  immediately  upstream  of  the  blade 
row  to  the  change  in  vorticity  across  the  blade  row.  It  is  expressed  in 


2»* 


finite  difference  form  Toy  approximating  the  time  derivatives  by  a  forward 
difference  scheme  and  the  spatial  terms  by  a  central  difference  apprcxima- 


and  again  using  Eq.  (38), 

(4 

which  must  be  satisfied  by  the  vorticity  at  the  blade  row  exit  plane. 

Flow  Field  at  Exit  Plane 

The  final  boundary  condition  states  that  the  static  pressure  is 
constant  along  the  flow  field  exit  plane,  CF.  From  Eq.  (23)  this  is  given 


COUPLED  PLOW  FIELD  AND  BLADE  ROW  ANALYSIS 


The  combination  of  flow  field  and  blade  passage  solutions  described  in 
the  previous  two  sections  was  programmed  for  the  UARL  UNIVAC  3108  digital 
computer  and  a  number  of  cases  were  examined.  In  this  portion  of  the  study 
the  system  geometry  and  the  loss  and  turning  characteristics  were  identical 
to  those  in  Ref.  7* 

A  cosine  wave  inlet  velocity  distortion  was  chosen  for  the  first  cases 
examined,  given  by  the  formula  (U-U*)/lf*  «  -0.3ccs6.  It  has  a  magnitude 
equal  to  30  percent  of  the  undistorted  velocity  (based  on  the  circumferential 
axial  velocity  average,  U*)  and  is  imposed  £  blade  row  circumferences  upstream 
of  the  blade  row.  This  is  shown  as  the  dashed  line  in  Fig.  10  for  an  inlet 
angle  of  42  deg.  Also  shown  in  this  figure  is  the  axial  velocity  dis¬ 

tribution  at  the  rotor  inlet  at  a  dimensionless  time  t*  *  2.94  after  the 
Imposition  of  the  distortion.  (Note  that  when  t*  s  Qt/2n  •  1  the  rotor  has 
made  1  revolution.)  The  same  steady-state  distribution  was  obtained  for  all 
values  of  t*  after  the  passage  of  the  initial  transient.  This  figure  shows 
that  the  distortion  pattern  is  attenuated  axially  and  shifted  circumferen¬ 
tially  in  the  direction  cf  rotation  by  the  presence  of  the  downstream  blade 
raw.  At  this  inlet  angle  the  rotor  blade  incidence  angle  is  below  the  stall 
angle,  5^  *  35  deg,  and  no  flow  instability  is  present.  (Note  that  the 
incidence  angle  can  be  defined  as  the  difference  between  the  blade  row 
stagger  angle  and  the  inlet  angle,  a  -  0*  -  01  .  In  this  case  0j*  =  63*5 
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deg  so  0^  -  42  deg  in  Fig.  9  corresponds  to  *  *  21.5  deg  and  the  stalling 
angle  0^  «  35  deg  is  equivalent  to  or  ■  28.5  deg.) 

The  effect  of  increasing  incidence  angle  is  shown  in  Figs.  11  through  13 
for  the  same  cosine  wave  distortion.  In  the  cases  examined  here  the  flow 
rate  was  held  constant  and  the  wheel  speed  was  increased.  This  has  the 
effect  of  decreasing  g1  and,  hence,  increasing  a  .  In  Fig.  11  a  2  deg 
increase  in  incidence  angle  causes  a  significant  flow  instability  to  occur 
in  the  low  velocity  region  (270  deg  <  0  <  360  deg),  as  evidenced  by  the 
lack  of  repeatability  for  the  three  values  of  t*  selected  for  this  figure. 

The  instability  appears  to  be  attenuated  as  the  high  velocity  region  is 
approached,  but  the  velocity  profiles  are  not  completely  repeatable  —  they 
merely  exhibit  similar  behavior  as  the  high  velocity  region  is  approached. 

In  Fig.  12,  for  a  «  25.5  deg,  the  unsteadiness  is  much  more  pronounced 
and  now  persists  over  a  larger  region  of  the  circumference.  It  should  be 
noted  that  the  fluctuation  in  velocity  at  the  blade  row  in  the  low  velocity 
region  is  now  of  the  same  order  of  magnitude  as  the  original  velocity  distor¬ 
tion.  Specifically,  note  the  change  in  the  nonnalized  axial  velocity  for 
t*  =»  2.42  and  t*  ■  2.90  at  0  »  270  deg.  As  before,  the  unsteadiness  appears 
to  attenuate  as  the  high  velocity  region  is  approached. 

Finally,  in  Fig.  13,  we  have  what  appears  to  be  a  completely  unsteady 
flow  which  persists  over  the  entire  circumference.  Waves  are  formed  at  the 
inlet  face  which  propagate  nonunifcrmly  in  the  direction  of  rotation  but  at 
a  fraction  (averaging  £)  of  rotor  speed,  indicative  of  rotating  stall.  This 
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is  shown  in  more  detail  in  Pig.  1^  which  i»  a  computer-generated  plot  for 
the  name  case  as  in  Pig.  13 •  The  horizontal  axis  is  again  circumferential 
position  and  the  vertical  axis  is  both  normalized  velocity  and  nondimens ional 
time.  The  first  curve  (at  the  bottom  of  the  figure),  denoted  "cosine  wave" 
iG  the  original  distortion  imposed  £  circumference  upstream  of  the  blade  row. 
The  remaining  curves  are  the  velocity  profiles  at  the  blade  row  inlet  taken 
at  a  series  of  values  of  t*  which  increases  from  bottom  to  top.  This  pseudo- 
three-dimensional  figure  shows  the  formation  of  a  velocity  deficiency  in  the 
low  velocity  region  which  propagates  in  the  direction  of  rotation  (to  the 
left)  and  attenuates  as  it  passes  through  the  high  velocity  region,  only  to 
reappear  and  strengthen  again  in  the  low  velocity  region. 

It  should  be  noted  that  in  Ref.  7,  which  analyzed  this  blade  row  geometry 
with  no  inlet  distortion,  rotating  stall  first  occurred  at  0^  =  35  deg  or 
o  =  28.5  deg,  which  is  a  higher  incidence  angle  than  any  shown  here.  With 
the  imposition  of  a  3°  percent  distortion  the  first  evidence  of  rotating 
stall  is  found  in  Fig.  12  for  a  »  25*5  deg  and,  to  some  extent  in  Fig.  11 
for  a  =  23-5  deg.  Thus,  the  presence  of  an  inlet  distortion  serves  to  reduce 
the  equivalent  steady-state  incidence  angle  at  which  rotating  stall  occurs, 
in  this  case  by  as  much  aB  5  degrees.  The  analysis  predicts  the  same  effect 
as  that  observed  in  rotor  experiments;  namely,  the  stable  operating  range  of 
the  rotor  is  reduced  by  the  presence  of  an  inlet  distortion. 

The  influence  of  the  unsteady  kinetic  energy  term  within  the  blade 
passage  is  examined  in  Figs.  15  and  16  by  setting  C^/2 nR^*  0.  This  has 


28 


the  effect  of  eliminating  the  entire  kinetic  energy  term  in  Eq.  (21). 

The  same  upstream  cosine  wave  distortion  was  used  as  in  Fig.  10.  A  com¬ 
parison  of  Figs.  10  and  15  for  a  =>  21.5  deg  and  Figs.  13  and  16  for  at  = 

27.5  deg  shows  that  the  absence  of  the  kinetic  energy  term  has  virtually  no 
effect  on  the  predicted  velocity  profile  at  "both  low  and  high  incidence 
angles.  This  is  in  apparent  contradiction  with  the  results  described 
earlier  in  the  section  entitled  Application  of  Blade  Row  Analysis,  in  which 
the  kinetic  energy  effect  was  dominant.  However,  in  the  earlier  results  the 
rotor  blade  row  was  preceded  by  inlet  guide  vanes  which  retarded  the  circum¬ 
ferential  attenuation  of  the  relatively  sharp  edged  distortion.  Hence,  the 

factor  dU/  dt  in  Eq.  (21)  was  important  in  this  application  of  the  analysis 

♦ 

and,  accordingly,  so  was  the  kinetic  energy.  In  the  present  case  the  cosine 
wave  has  a  much  smaller  circumferential  gradient  to  begin  with,  and  the  axial 
attenuation  further  reduces  the  influence  of  the  velocity  time  derivative, 
dU/dt. 

Next,  a  6teeper  and  narrower  distortion  pattern  was  examined,  shown  by 
the  dashed  line  in  Fig.  17 •  In  this  case  the  wheel  speed,  Q  ,  and  the  average 
axial  velocity,  U*,  were  the  same  as  in  Fig.  10.  This  time  the  distortion 
cosine  wave  extended  only  over  the  range  90  deg  <  8  <  270  deg,  and  was 
flat  elsewhere  (l.e.,  this  was  a  half-wave  180  deg  distortion).  This  is 


given  by 


0.15,  O<0<9Odeg,27Odeg  <6  <  360deg 

(43) 

0.3cos2|(9-  -Jj  -0.15 ,9Odeg<0<  270deg 

This  waveform  and  the  full  circumference  wave  used  earlier  both  have  the 
6ame  circumferentially  averaged  mean  value  of  axial  velocity  and  both  have 
the  same  peak-to-peak  amplitude. 

It  is  readily  seen  that  this  distortion,  extending  over  only  180  deg 
of  the  circumference,  has  a  higher  harmonic  frequency  content  than  the 
previous  distortion  which  extended  over  the  full  3&0  deg  of  circumference. 

In  Fig.  17,  for  t*  =  2.94,  the  peak-to-peak  distortion  at  the  blade  row 
inlet  (solid  curve)  is  smaller  In  magnitude  than  in  Fig.  10.  It  is  believed 
that  this  amplitude  attenuation  takes  place  because  the  inertia  of  the  flow 
field  prevents  a  rapid  response  to  the  higher  harmonics.  In  effect,  the 
fluid  inertia  behaves  as  a  low  paSB  filter,  permitting  the  response  to 
follow  the  lower  harmonics  of  the  upstream  distortion  but  suppressing  the 
response  at  the  higher  harmonics.  As  before,  the  distortion  at  the  blade 
row  inlet  is  shifted  in  the  direction  of  rotation. 

Figures  18  and  19  show  the  effect  of  increasing  incidence.  Again  the 
flow  instability  or  rotating  stall  develops  in  the  low  velocity  region  and 
attenuates  in  the  high  velocity  vegion,  with  the  effect  becoming  more 
pronounced  as  a  is  Increased .  However,  this  time  the  instability  is  unable 
to  grow  as  before  around  the  entire  circumference,  even  at  the  highest 


incidence  angle,  probably  because  of  the  inertia  effects  cited  above. 

Finally,  in  Fig.  20,  the  effect  of  suppressing  the  kinetic  energy 
term  is  again  examined  for  this  higher  relative  frequency  by  setting 
Cax/2nRM  *  0.  The  upper  curve  is  for  low  incidence,  a  *  21.5  deg,  and 
the  lower  curve  is  for  high  incidence,  a  «  27*5  deg.  A  comparison  of 
the  upper  part  of  Fig.  20  with  Fig.  17  for  a  ■  21.5  deg,  and  the  lower 
part  of  Fig.  20  with  Fig.  19  for  a  *  27*5  deg  again  shows  a  minimal  in¬ 
fluence  of  the  kinetic  energy  term  on  velocity  response.  Here  it  is 
believed  that,  although  the  initial  distortion  is  steeper  than  before,  the 
axial  attenuation  that  takes  place  at  the  blade  row  inlet  reduces  the 
contribution  of  the  factor  3U/  at  to  approximately  the  same  level  as  in 
the  full  circumference  distortion. 


RESULTS  AND  CONCLUSIONS 


The  results  of  this  investigation  are  as  follows: 

1)  A  nonlinear,  large  disturbance  theory  has  been  developed  which 
couples,  interactively,  the  flow  through  the  blade  passages  of  a  turbo- 
r.achine  blade  row  and  the  flow  through  the  surrounding  field. 

2)  The  theory  predicts  the  degradation  of  the  stable  operating  range 
of  an  isolated  rotor  due  to  the  presence  of  an  upstream  distortion. 

3)  Use  of  the  uncoupled  blade  passage  theory  shows  that  for  a  sharp- 
edged,  unattenuated  distortion  the  kinetic  energy  term  is  dominant. 

4)  Use  of  the  coupled  theory  and  cosine  wave  distortions  Bhcws  that 
the  kinetic  energy  term  is  less  important  because  of  the  reduced  circum¬ 
ferential  gradient  of  velocity  and  that  the  dominant  effects  arise  from  the 
inclusion  of  total  pressure  loss  and  from  the  coupling  between  blade  row 
pressure  changes  and  flow  field  vorticity  changes.  However,  if  sharp  velocity 
gradients  exist,  the  kinetic  energy  term  would  probably  contribute  strongly 

in  these  regions  of  high  gradient. 

The  conclusions  reached  in  the  course  of  this  study  are  as  follows: 

1)  The  predicted  degradation  of  the  stable  operating  range  of  a  single 
blade  row  due  to  the  presence  of  an  axial  distortion  is  of  the  same  order  as 
the  available  measured  results  from  multistage  machines . 

2)  At  low  incidence  angles  an  upstream  axial  distortion  is  attenuated 
spatially  as  it  approaches  a  blade  row,  and  the  flow  through  the  blade  row 


remains  stable. 

3)  At  high  incidence  angles  an  upstream  axial  distortion  having  a 
circumferential  wave  length  equal  to  the  annular  circumference  will  produce 
d isturbances  in  the  region  of  the  blade  row  which  are  amplified  axially  in 
the  low  velocity  region  and  attenuated  in  the  high  velocity  region. 

4)  These  disturbances  propagate  circumferentially  as  a  rotating  stall 
in  the  direction  of  rotation  at  approximately  one-half  the  wheel  speed  of 
the  rotor. 

5)  An  upstream  axial  distortion  which  occupies  one-half  the  annular 
circumference  will  be  attenuated  more  than  the  longer  wave  length  distortion 
because  the  inertia  of  the  fluid  restricts  the  response  capabilities  of  the 
fluid  system  to  higher  harmonics. 

6)  For  the  one-half  circumference  distortion  at  high  incidence  angles 
the  tendency  for  the  creation  of  rotating  stall  cells  is  reduced. 


aMagggas. 
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LIST  OF  SYMBOLS 


Area,  ft2 

Axial  component  of  "blade  chord,  ft 

Functional  form  of  vorticity,  (ft/see)/ft 

Functional  form  of  stream  function,  ft2/sec 

Unit  normal  vector 

Static  pressure,  lb/ft2 

Total  pressure,  lb/ft2 

Total  pressure  loss,  lb/ ft2 

Resultant  velocity,  ft/sec 

Flow  rate  through  blade  passage,  ft3/sec 

Axial  flow  rate  through  entire  flow  field,  ft^/sec 

Mean  radius  of  blade  row  annulus,  ft 

Time,  sec 

Dimensionless  time,  *  Qt/2rT 
Axial  velocity,  ft/sec 

Circumferential  average  of  Axial  velocity,  ft/sec 
Volume,  ft3 

Tangental  velocity,  ft/sec 

Axial  coordinate,  ft 

Tangential  coordinate,  ft 

Cascade  Incidence  angle,  »  8£  -  ,  deg 
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LIST  OF  SYMBOLB  (Continued) 


Dimensionless  time  constant 


01,  02 


Inlet  and  exit  flow  angles,  aeg 


5  metal 


Cascade  stagger  angle,  deg 
Vorticity,  (ft/sec)/ft 
Circumferential  coordinate,  rad  or  deg 
Local  flow  angle  within  "blade  passage,  rad 
Local  angle  of  "blade  camber  line,  rad 
Chordwise  average  of  blade  camber  line  angle 
Fluid  density,  lb  sec^/ft^ 

Stall  angle  parameter,  ratio  of  instantaneous  incidence 

angle  to  stall  angle 

Tangential  gap,  ft 

Unsteady  loss  parameter 

Steady-state  loss  parameter 

Stream  function,  ft^/sec 

Rotor  wheel  speed,  rad/sec 


Special  symbols,  superscripts  and  subscripts 


Iaplacian  operator 
Gap  average 


Blade  row  inlet  plane 


Blade  row  exit  plane 


37 


CIRCUMFERENTIAL  POSITION, 


CIRCUMFERENTIAL  STATIC  PRESSURE  PROFILE  AT  OUTER  DIAMETER  FOR 
180  DEG  DISTORTION  EXTENT  -  PEAK  PRESSURE  CONDITION 


CIRCUMFERENTIAL  POSITION 


CIRCUMFERENTIAL  STATIC  PRESSURE  PROFILE  AT  OUTER  DIAMETER  FOR 
180  DEG  DISTORTION  EXTENT  -  NEAR  SURGE 


CIRCUMFERENTIAL  POSITION, 


CIRCUMFERENTIAL  POSITION 


SCHEMATIC  DIAGRAM  OF  COMPRESSOR  FLOW  FIELD  WITH  SEMI-ACTUATOR  DISK  REPRESENTATION  OF  BLADE  ROW 


NORMALIZED  AXIAL  VELOCITY  COMPONENT, 


FIG.  10 


EFFECT  OF  COUPLING  ON  360  DEG  COSINE  WAVE  DISTORTION  INCLUDING  KINETIC  ENERGY 


TERM  FOR  fa  =  42  DEG 
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NORMALIZED  AXIAL  VELOCITY  COMPONENT,  'U-U*)/U 
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EFFECT  OF  COUPLING  ON  360  DEG  COSINE  WAVE  DISTORTION  INCLUDING  KINETIC  ENERGY 
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NORMALIZED  AXIAL  VELOCITY  COMPONENT.  (U-UM/U 
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EFFECT  OF  COUPLING  ON  180  DEG  COSINE  WAVE  DISTORTION  WITH  KINETIC  ENERGY  TERM 
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